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Abstract 

o , 

I A lattice quantum gravity model in 4 dimensional Riemannian spacetime 

is constructed based on the SU(2) Ashtekar formulation of general relativity. 
This model can be understood as one of the family of models sometimes called 
"spin foam models." A version of the action of general relativity in continuum 
is introduced and its lattice version is defined. A dimensionless "(inverse) cou- 
pling" constant is defined so that the value of the action of the model is finite 
per lattice point. The path integral of the model is expanded in the characters 
and shown to be written as a sum over surface-like excitations in spacetime. A 
■ 3 dimensional version of the model exists and can be reduced to lattice BF the- 

ory. The expectation values of some quantities are computed in 3 dimensions 
q | and the meanings of the results are discussed. Although the model is studied 

\Q • on a hyper cubic lattice for simplicity it can be generalized to a randomly 

triangulated lattice with small modifications. 
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1 Introduction 



A surface theoretic view of non-perturbative quantum gravity was studied some time 
ago based on 3 dimensional general relativity as a reformulation of the Ponzano-Regge 
model in Riemannian spacetime [I| . Since 3 dimensional general relativity has no local 
degrees of freedom, the Ponzano-Regge model is known to be a topological field theory 
called lattice BF theory in 3 dimensions. A 4 dimensional generalization was made 
by Ooguri and is known as lattice BF theory in 4 dimensions, also a topological field 
theory 0. This Ooguri model also allows the surface theoretic view. However, since 
4 dimensional general relativity contains local degrees of freedom, the Ooguri model 
is not a quantization of general relativity. 

In order to pursue 4 dimensional generalizations of the surface theoretic view of 
the Ponzano-Regge model with local degrees of freedom, some investigations have 
been being made. They can be divided roughly to four strategies. One [||] is the use 
of canonical theory called loop quantum gravity |4[] based on Ashtekar formulation 
of general relativity ||. It uses a Hamiltonian of the canonical theory to compute 
quantum amplitudes. In this strategy, the definition of "correct" Hamiltonian is still 
uncertain although the canonical theory is able to provide physical interpretations. 
Another || [7], [8| is the use of the Plebanski action of general relativity in Reiman- 
nian spacetime. This action has the form of BF theory with additional constraint 
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terms M. It computes a path integral of the action. In this strategy, a well regu- 
lated computation of the constraint terms is still missing and, in addition, the action 



functional is available only in Reimannian spacetime. Another |L0[ [1 1] is the use of 
"quantized" lattice BF theory and additional "quantum level" constraint conditions. 
In this strategy, the formulation is strictly mathematical and it necessitates help from 
other theories for physical interpretations. The other 0] is based on physical rea- 
soning such as causality and statistical criticality. This strategy wants a systematic 
formulation. 

In the present work, a fifth strategy is applied to construct a surface theoretic 
model for 4 dimensional quantum gravity. We introduce a version of the action of 
general relativity in continuum. This action is derived from the action proposed some 



time ago by Samuel and by Jacobson and Smolin The Ashtekar formulation 

can be constructed from the latter. We define a lattice version of this action and 
its path integral. The variables integrated are closely related to the SU(2) (real 
valued) variables of the Ashtekar formulation. However, the model avoids the use of 
Hamiltonian of canonical theory, which is a main difficulty so far for quantization. The 
spacetime the model assumes has signature (+1, +1, +1, +1), Riemannian spacetime. 
In order to make the lattice action finite per lattice point, a dimensionless "(inverse) 
coupling" constant is defined in analogy to Wilson's formulation of lattice gauge 
theory []14" |. The finiteness of the action with finite lattice allows the path integral 
without gauge fixing finite. The path integral of the model is computed as a character 
expansion and shown to be written as a sum over surface-like excitations in spacetime. 

We show that a 3 dimensional version of the model exists and can be reduced to 
lattice BF theory. In this sense, the model is a 4 dimensional generalization of the 
Ponzano-Regge model with local degrees of freedom. The expectation values of four 
quantities are computed in 3 dimensions. Two of them are basic variables of the theory 
and are SU(2) gauge dependent quantities. Another is the Wilson loop. It is SU(2) 
gauge independent but diffeomorphism dependent. The other is the action itself. It is 
invariant under not only SU(2) gauge transformations but also diffeomorphisms. The 
first three are shown to be always zero unless an appropriate gauge is fixed. The other 
is zero if the (inverse) coupling constant is brought to an appropriate limit. These 
are very simple results but important checks for consistency. We discuss the meaning 
of these results. Note that this kind of computations have not been performed in 
the other strategies while the computation of the path integral (or sometimes called 
transition amplitude, projector or partition function) is always focussed on in all the 
strategies. 

The model is defined on a fixed hyper cubic lattice for simplicity. However, the 
construction can be generalized to a randomly triangulated lattice with small modifi- 
cations. Further more, the sum over different triangulations can be done, if necessary, 
by applying a technology developed in a recent work [15]. 

This paper is organized as follows. In Sec. & we restructure the action functional 
of general relativity so that it can be utilized to define the model. In Sec. the model 
is defined and its path integral is studied. In Sec. |4], three dimensional version of the 
model is studied. In Sec. [5], we conclude the work. In Appendix |A], mathematical 
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formulae utilized in this work are listed. 



2 The action of general relativity 

The action we make use of is the one proposed by Samuel and by Jacobson and 
Smolin fl3| . It is the Palatini action but an additional term. The additional term 
is a topological term and not responsible to the local degrees of freedom of general 
relativity. From the Palatini action, one can construct a canonical theory by splitting 
space and time. However, it contains constraints of the second class in addition to 
those of the first class in the sense of Dirac. Therefore, the constrained phase space 
is not well defined unless the second class constraints are explicitly solved. 

On the other hand, from the present action, one can construct Ashtekar's canonical 
theory, which contains only constrains of the first class. It means that the constrained 
phase space is well defined in the sense that it is invariant under transformations 
produced by Hamiltonian and diffeomorphism generators. The difference between 
the two canonical theories is due to the canonical transformation produced by the 
topological term in the Samuel- Jacobson-Smolin action. We restructure this action 
for the use in the following section. 

Let us introduce variables we need. We use tetrad field el rather than metric 



7< 



g^ u . They are related with each other by = e e^rjij. Here rju is the Euclidean 
or Minkowskian metric of "internal" spacetime depending on the value of a with sig- 
nature (a, +1, +1, +1). The capital alphabets I, J ■ ■ ■ are used for internal spacetime 
indices {0, 1, 2, 3} and the Greek letters are for spacetime indices {0, 1, 2, 3}. Below, 
we also use the lower case alphabets i,j--- for the space componets of the internal 
spacetime indices {1,2,3} and for the adjoint indices of SU(2) group. 

The so-called spin-connection u IJ is written in terms of the tetrad field as follows. 



Here e^ 1 is the inverse tetrad. The Riemann tensor R,r„ is written in terms of the 
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spin-connection as follows. 

Rjl := 2d[Xf + HV" ( 2 ) 

In terms of these traditional variables, we define connection and curvature variables 
Aj/ and F 1 ^ as follows. 

K J := ^ + \^ KL ^ = \^ KL A™ + (1 - AX (3) 



2' e KL^iiv — 2 

Kt = 2 <9[X] + AffA v]K J + (1 - T^V]* ■ « - v » 



Kl ■= R % + o^ U klR^ L = -^ 1J klF« L + (1 - 7 2 *)i#, (4) 
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Here 7 is a parameter and chosen to be 7 = 1 to relate the model with SU(2) (real 
valued) Ashtekar formulation. If one chooses the value such that 7V = 1, then the 
variables Aj/ and are self-dual with respect to the internal spacetime indices. 
We keep 7 unspecified in the course of development within this section. 

Among the components of Aj/ and F™, we use specifically A / and is 
written in terms of A®* and tuj/ as follows. 

F% = 2d^A§ - a^^Af + (1 



7V) 



O, ,0fc, , % I "^,_Ufc , ,mn, , i 
ZuJ [^v}k + q7 6 mn W \p, ^ u ]k 



.(6) 



The pull-back of A / to foliated space of spacetime is precisely the SU(2) (real valued) 
connection variables of the Ashtekar canonical formulation. 

Let us explicitly write the action. First, the Palatini action derived from the 
Hilbert-Einstein action is 



S [e,u] := 4 / d 4 Xy/agR = — f d 4 xee^ejR 
l p J l p J 

= wl d4 ^ a tTJK L ey u Rl L . 



u 

/J.U 



(7) 



Here, g is the determinant of metric g^, R is the scalar curvature, l p is the Planck 
length constant, e is the determinant of ej and e^ uXa is the alternating density. By 
adding a topological term as mentioned above, the action we make use of is 



S[e,A] :=-Jd 

= _L / d'xe^e^jF^. 
7% J 



7f1 e KL K \ 



KL 
a 



R 



(8) 



This is the Samuel- Jacobson-Smolin action. By dividing space and time indices for 
internal spacetime and writing in terms of e°, and A ^, the action becomes 



S[e,A] = — / d 4 xe^ [2e,„e„^ + e^F^ 
7'p J 

= ^2 / d ' xe ^ [^e m Fl + e, t e UJ [ 7 e« Qk F° x k a + (1 - 7 V)J2& 
= ^| / [* + 7 eVi e *) ^ + (1 - T^e^-H 



+(i- 7 V) 



a (2^ + 7 eVl e *) (2o;ga; ff]Jfc * + -76 0fe „ 



171 Ti 



(9) 



Here, we consider uf/ and if^ as functions of e° u and e„. Note that if 7V = 1, then 



the terms proportional to (1 — 7 2 c) go away and the action gets simplified. We will 
utilize the case that 7 = 1 and a — 1, Riemannian spacetime. 
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3 The model 



3.1 The action 

The action discussed in the previous section for Riemannian spacetime (a = +1) with 
7 = 1 is 

S*[e, A] = ±f d'xe^ + e\ jk e^ (2^4] " ^ jt A4) ■ 

Hereafter, we write A 1 ^ for We note that if 2e^eJ,] + e^e^ i s replaced by 

a 2-form variable denoted by with the condition that & vXa B l B Xa is position- 
wisely proportional to 5 lJ , then the resulting action is precisely the Plebanski action. 
This additional condition is supposed to restore B^ u = 2ej^e^ + e 0t '^e^e*. Ways of 
constructing a surface theoretic model starting from the Plebanski action have been 
explored. The condition is usually imposed with Lagrange multipliers. However, it 
has turned out that the inclusion of the condition in terms of Lagrange multilpiers is 
still a difficult task. In the present work, we do not follow any of those ways. Instead, 
we define a corresponding action directly on the lattice without introducing B l and 
utilizing Lagrange multipliers. 

In order to define the model, we rewrite the action in a slightly compact fash- 
ion so that the construction of the lattice version is straightforward. An idea is the 
following. The tetrad variables (or their modifications) are often considered su(2) 
algebra valued. For example, in the Ashtekar formulation, the pull-back of the (in- 
verse) tetrads to space have su(2) algebra indices and , in the Plebanski formulation, 
the 2-form variable B % are su(2) algebra valued. However, we do not consider here 
the tetrad variables su(2) algebra valued but do consider them SU(2) group valued. 
If one defines e M := e° + icr^e 1 with the Pauli matrices <jj (i — 1, 2, 3), then the com- 
plicated expression 2e^e l v ^ +€ 0l jk e3 j, e t becomes just the non-trace part of e^e v . If e M is 
normalized so that e°e° + e*e* + e^e^ + e^e^ = 1, then it is an SU(2) group element. 
This small change of view dramatically simplifies the construction of the model and 
makes the extension to other dimensional spacetimes straightforward. 

Before introducing the lattice, let us count the number of degrees of freedom of 
the variables, and A 1 ^ have respectively 16 and 12 degrees of freedom per spacetime 
point. From thses variables, if one constructs canonical pairs of variables, one finds 

2 physical degrees of freedom from and 2 from A 1 . The other degrees of freedoms 
are for the constraints and gauge trajectries in addition to the freedom of specifying a 
spacetime foliation. Let us be slightly more concrete. Out of 16 degrees of freedom of 
ejt, 3 fix a spacetime foliation and 4 are used for Lagrange multipliers imposing three 
diffeomorphism and a Hamiltonian constraints. Out of 12 degrees of freedom of A 1 , 

3 are used for Lagrange multipliers imposing three Gauss gauge constraints. Hence, 
9 degrees of freedom in each are left. These are exactly the number of degrees of 
freedom of the Ashtekar canonical theory. Then by subtracting 7 for the constraints 
and 7 for the gauge degrees of freedom, one finds the correct number of physical 
degrees of freedom. 
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3.2 Lattice 



We fix a lattice. We do not take sum over different lattices. In general, the model can 
be defined on a randomly triangulated lattice. However, technically it is difficult, if 
not impossible, to compute the quantum amplitudes analytically or even numerically 
if the structure of lattice varies from place to place. For this reason, we use a hyper 
cubic lattice. First we introduce a pair of hyper cubic lattices. One is dual to the 
other. Denote them A and A* respectively. Then we further restrict ourselves to 
the case that A* is A itself. This is possible because the dual lattice of a hyper 
cubic lattice is a hyper cubic lattice. This restriction dramatically simplifies the 
construction of the model and helps computations in practice. 

The lattices A and A* are a pair of 4-dimensional hyper cubic lattices dual to 
each other. They consist of 0-cells (vertices), 1-cells (edges), 2-cells (square faces), 
3-cells (cubes) and 4-cells (hyper cubes). Every k-cell of one lattice intersects with a 
(4-k)-cell of the other lattice at a point inside the cells. Every cell of the lattices has 
one and only one intersection. The correspondence of k-cells of A and (4-k)-cells of 
A* is one-to-one. The corresponding k-cell of A and (4-k)-cell of A* are called dual 
to each other. There is a limit at which every intersection coincides with one of the 
vertices so that the two lattices are identical A = A*. 

Let us be more specific. Fix A and A*. Both have a lattice spacing e. One is the 
displacement of the other by the distance e/2 in all the four directions. Let x denote 
(the position of) a vertex of A. The (positions of) adjacent vertices are x ± with 
fx — 0,1, 2, 3, where is the increment in the /z-th direction with the magnitude e. 
Note that x ± e°, for example, is a short hand notation for (x° ± e, x 1 , x 2 , x 3 ). Let x* 
denote (the position) of a vertex of A* at x + e°/2 + e 1 /2 + e 2 /2 + e 3 /2. The face of 
A specified by (the positions of) four vertices x, x + , x + e u and x + e M + e v is dual 
to the face of A* specified by (the positions of) four vertices a and 

x* — e x — e a , where the indices are chosen such that e^ vXa = 1. There are six pairs of 
faces for given x up to the orientation of the face. These dual faces are important in 
our construction of the model. After taking the limit A* — > A such that x* coincides 
with x, the dual faces are still clearly defined on the single lattice A. The intersection 
of a dual pair of faces, one basing at x and the other at x* , coincides with x = x* 
at the limit. Notice that the face basing at x* and in the fiu plane does not coincide 
with the face basing at x and in the same plane at the limit. One is in the /i and v 
directions while the other is in the —\i and —v directions from the coincident base 
point x = x*. We will define our model on the lattice A(= A*). 

3.3 The lattice action 

We define an action on the lattice A so that it converges to the action in continuum 
as the lattice spacing goes to zero. First, let us define variables on the lattice as 
follows. For an edge basing at x and pointing at the /z-th direction, define 

Ux) := ^(x)+ia i q i (x) := eM^(x)a t /2], (11) 
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rj^x) := rf^x) + iatf^x) := {^%)- l e[e%x) + ie^x)^]/ p^(x), (12) 
p,(x) := (^Xr'e^xW + \el(x)\z + \el(x)\* + \<*(x)\*. (13) 

Here, e is the lattice spacing and cr^'s are the Pauli matrices, that is, 





'0 


r 




'0 -% 




'1 " 


<7i := 


.1 


o_ 


, ^2 := 


J o _ 


, ^3 : = 


.0 -1. 



Note that the inside the square root in p M is equal to the metric element g^. (3 is the 
upper bound of l~ 2 s 2 g^(p = 0, 1, 2, 3) and introduced to regularize the path integral 
defined below. p M is no n- negative and less than 1, and and r]^ are SU(2) matrices 
normalized such that |C°| 2 + |C*| 2 + \( 2 \ 2 + IQ 2 = 1 and |r?°| 2 + |^| 2 + |?? 2 | 2 + \vfc\ 2 = 1- 
For a face basing at x and enclosed by four edges in the p and z/-th direstions, define 

U^{x) := Ux)Ux + e^dix + e")Cj(x), (15) 

without summing over p and v. Here, e M is the increment in the p-th direction with 
the magnitude e. 

The action on the lattice A is then defined as follows. 

Sa[(,V,p\ := -P E ^ x ° P ,(x) Pu (x)Tt[ V 1(x) Vu (x)U x ^)]- (16) 

Here, the sums over p, u, A and a have been performed. This lattice action is finite 
on the lattice A with e and (3 fixed if the lattice size is finite. This is because and 
r]^ are compact SU(2) variables and p^ has lower and upper bounds in addition to 
the fact that the number of degrees of freedom is finite on the lattice if the lattice 
size is finite. The finiteness of the action lets the path integral without gauge fixing 
non-divergent. In other words, the path integral along a gauge trajectory produces 
just a finite overall multiplicative constant. We call (3 "inverse coupling" in analogy 
to Wilson's formulation of lattice gauge theory. 

Let us examine if the lattice action really converges to the continuum action as 
the lattice spacing e goes to zero. Take the limit e — > 0. In the limiting process, it 
can be shown that 

Pp^x)p v {x)r]l(x)r] v] (x) = l; 2 e 2 ia, t (^e^e^x) + e 0i jk ^{x)e h v {x)) , (17) 
U [Xa] (x) - (2d [x Ai } (x) - e 0i jk A{(x)A k a (x)) + G(e 3 ). (18) 

Therefore, the lattice action converges to 

Sa[(, v, p] - K 2 E ^ vXa (K( X M X ) + e V J >)*)) x 

(2d x Al(x) - e 0i jk A{(x)A k a (x)) + 0(e 5 ). (19) 
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3.4 The lattice path integral 

The exponential of the action can be considered as (an extension of) the graph- 
cylindrical function discussed in [[L6[]. The integral measure we use for the variable 



A 1 ^ is the Ashtekar-Lewandowski measure [16]. The measure for SU(2)-compactified 
part of e° and is analogously defined and that for the other part of them is the 
one discussed in Pfl. The path integral has the form of 



dp(A)dp(e)^ Aii ,(A, e) ■= / II H d C^( x )dVn{ x ) d Pt( x ) x 

J zGA At 

VKCOdj A), • • • C0«, 4); 77(3:1, e), • • • r/(x m , e); p(x h e), • • • p(x n , e)), (20) 

where ^a,^ is a cylindrical function defined on the lattice A in terms of a complex 
valued integrable function ip on [577(2)]' x [SU(2)] m x [0, l] n . d(^ and are the 
Haar measures for SU(2) elements and 77^ respectively dp* = Adp^p 3 ^ is the radial 
integral of the 4-dimensional unit sphere. p M is integrated from to 1. 

The path integral we will compute in the following section is defined as follows. 

Z A := J c/C^^p 4 e" i5A[C '' ? ' p] 

= f d(dr]dp 4 TT eifc^p^pA^lvtWvAxWxAx)]^ ^1) 

Note that we integrate the exponential oscillation form e lS rather than the expo- 
nential decay form e~ s . The latter is commonly and successfully used for Euclidean 
quantum field theories whose action has quadratic form. However, the action of gen- 
eral relativity is not quadratic but linear in each of the variables. Because of this fact, 
the action is not bound from the below and the use of the exponential decay form 
runs into serious technical problems. In addition, it unlikely contains the classical 
limit since the value of the action for classical solutions is zero. The former is used in 
other strategies in quantum gravity [[[], 0, |3], || [7|, |^, [LT], |TT], |TJ] even for Riemannian 
spacetime. In the exponential oscillation form, large magnitudes of the action less 
contribute to the path integral. In order to examine possible consequences of the use 
of the exponential decay form, we simply replace (3 by i/3 in the results of the use of 
the exponential oscillation form. 

In terms of the path integral, define the expectation value of an observable X as 
follows. 



POa := ^a 1 / dC<Mp 4 X[C, 77, p}e- iS ^>™\ (22) 



In particular, we are interested in four quantities: (77^), (Cn), (Trt/^) and (5a). The 
first two are basic variables of the theory and are SU(2) dependent quantities. They 
should vanish. This is because there is a symmetric structure along the SU(2) gauge 
trajectories. It is known that the expactation value of the basic variables of lattice 
QCD vanishes unless an appropriate gauge is fixed. This fact is often utilized to 
study gauge fixing ambiguities on lattice. The same thing should happen here since 
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the Gauss gauge constraint structure of the present theory is common with the SU(2) 
version of QCD. 

The third is (a smallest of) the Wilson loop. It is SU(2) gauge independent and a 
good physical observable for Gauss gauge invariant theories such as lattice QCD. Here, 
it is possible that the expectation value of the Wilson loop vanishs unless an appro- 
priate gauge is fixed. This would mean that the Gauss gauge invariance is not enough 
to be a physical observable for quantum gravity and a possible symmetric structure 
along diffeomorphism gauge trajectories annihilates the expectation value of the Wil- 
son loop. Quantum gravity is invariant not only under Gauss gauge transformations 
but also under spacetime diffeomorphisms. The Wilson loops are not diffeomorphism 
invariant and hence cannot be physical observables of quantum gravity. 

The last one is the action of the theory. It is invariant under SU(2) gauge transfor- 
mations and diffeomorphisms and the only known local physical observable in general 
relativity. It should vanish if the inverse coupling constant is brought to infinity but 
diverge in the same limit if the path integral has the exponential decay form. The 
former would be an indication that the path integral successfully eliminates the de- 
grees of freedom constrained since the value of the action on the constraint surface 
must be identically zero. The latter would be an indication that the path integral of 
exponential decay form cannot capture the classical solutions. These are non-trivial 
checks for consistency. These facts are examined in 3 dimensions. 

3.5 Character expansion 

We compute the path integral as a character expansion and show that it can be 
written as a sum over surface-like excitations in spacetime. First we rewrite the path 
integral as follows. 

Z A := J dC^P 4 e" i5A[C '"' pl 

= J d(dridp 4 JJ JJ e i2 f 3 Pv( x )p»( x ) Tr lvU x )vv(x)u Xa (x)] e -i2i3p^(x)p„ 

(23) 

Here, (puXa) denotes the even permutations of (0123); there are six terms, that is, 
(0123), (0231), (0312), (1203), (3102), and (2301). Then use the formula © to 
expand to the characters as follows. 

with 

(25) 



9 



Here, Xj{U) is the character of SU(2) element U in the spin-j representation and 
and l vft are spins taking values 0, ~, 1, | • • •. J m (x) is the Bessel function of the first 
kind. One of its definitions is given by 



From (|24|) , one can understand how surface-like excitations emerge. U\ a (x) is a 
parallel transport along the four edges enclosing a face of the lattice. The face is in 
the Acr-plane and bases at x. Let f\ a (x) denote this face. (We do not distinguish the 
orientations; f\ a {x) = f a \(x).) Associate the spin of Xl»v{v\{ x )Vv{ x )U\a{x)) to the 
face f\ a (x). In the same way, l UfM is also associated to the same face. By integrating 
r/^ at x, one introduces an intertwiner to combine six spins £ M „ and l vll [y ^ // for 
given /i) associated with f\ a (x) (A and a are such that ^ uXa = 1). By integrating 
Cfj, shared by six faces f^x) and f^ u (x — e u ) [y ^ [i for given /i), one introduces an 
intertwiner to combine the twelve spins associated with the six faces. Therefore, after 
integrating all 77's and £'s, the spins and the intertwiners combining the spins are left. 
Two spins are associated with a single face and two intertwiners are associated with 
a single edge. These facts are more clearly understood in the 3 dimensional version 
of the model as examined in Sec. The path integral is now written as a sum over 
spins and intertwiners associated with faces and edges respectively. A set of faces 
jointed by edges represent a surface-like object, mathematically called 2-dimensional 



piece- wise linear cell complex [fn]| . A surface-like excitation is a 2-dimensional piece- 
wise linear cell complex with faces labeled by non-zero spins and with edges labeled 
by intertwiners. A face with spin-0 means the absence of the face in a surface-like 
excitation. 

The integrations of 77's and £'s can be explicitly done by using the formulae (58), 
(j39p and fl40|). The resulting coefficiens contain the spins associated with faces and 
other spins representing intertwiners associated with edges. The coefficients are ex- 
pressed in terms of the so-called 3j-coefficients (pLt|). The integrations of p's are 
difficult because of the presence of two p's in the argument of the Bessel functions. 
This difficulty is absent in the 3 dimensional version of the model. 



4 The model in 3-dimensions 
4.1 The action and path integral 

From the action in 4 dimensions, it is straightforward to write its 3 dimensional ver- 
sion. We simply replace l~ 2 d i xe aX ^ l/ e^ex in the 4 dimensional action by /~ 1 <i 3 xe AAtI/ eA. 
We compute the corresponding path integral on the lattice and show that it is reduced 
to lattice BF theory in 3 dimensions in order to claim that 3 dimensional version of 
the model exists and is a correct quantum gravity model. 

The action of general relativity in 3 dimensional Riemannian spacetime is defined 
as follows. 

S R3 := y f d 3 xe x ^e x > - e^^X) • (26) 

try J 
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Here, the Greek letter indices are {0, 1, 2} for spacetime and lower alphabet indices 
are {1,2,3} for internal spacetime. 

Let us count the number of degrees of freedom. Out of 9 degrees of freedom of 
at each spacetime point, 3 are for Lagrange multipliers imposing two diffeomorphism 
and a Hamiltonian constraints. Out of 9 degrees of freedom of A 1 at each spacetime 
point, 3 are for Lagrange multipliers imposing 3 Gauss constraints. Hence, 6 degrees 
of freedom in each are left. They are cancelled by the 6 constraints and 6 gauge 
degrees of freedom and one finds no local degrees of freedom consistently with general 
relativity. Note that in canonical theory one sees the vectors specifying a spacetime 
foliation multiplicative e* . Therefore, the integral over all foliations is taken into 
account by the integrals of e* . 

The lattice action is defined as follows. 

S A3 [t lV} p] := -p £ e x ^p x (x)Tr[ Vx (x)U^(x)]. (27) 

Here, A3 denotes the 3 dimensional cubic lattice. Note that the definitions of the 
variables are the same as for the 4 dimensional model, (|TT|) , (|T2"D and (|T3"D, and the 
additional degrees of freedom rf^ cancel in the lattice action and do not play any role. 
Then the lattice path integral is defined as follows. 

Z A:j := [ dCdridp* J] e^^^^-^l 

= f d(dr]dp 4 J] n e ^^A^] e -^A^AtM (28) 

1EA3 (A M ")=(012), 
(120), (201) 

Then use the formula (|36|) to expand to the characters as follows. 



Za 3 = d(dr,dp 4 II II E ^x(fidhp})Xl^{VxU^)xhM U M 29 ) 

X£A 3 (A t »")=(012), la v ,l v a 
(120), (201) 

with 

UP,{1,P}) ■= {2l " V + 1)(2 V 2 1)( " 1)W " W ^ + i(2/?pA)J^ +1 (2 j 9 PA ).(30) 

P P\ 

Here, one can understand the emergence of surface-like excitations, in other words, 
2 dimensional piece-wise linear cell complexs with faces labeled by non-zero spins 
and with edges labeled by intertwiners. is a parallel transport along the four 
edges enclosing a face of the lattice. The face is in the /x^-plane and bases at x. Let 
ffj,is(x) denote this face. Associate the spin of Xi^iVxU^) to the face /^(x). In 
the same way, l Vil is also associated to the same face. By integrating i]\ at x, one 
introduces an intertwiner to combine two spins 1^ and l vll associated to /^(x) (p 
and v are such that e x ^ v = 1). By integrating shared by four faces f^x) and 
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ffiu{x — e") (v ^ /i for given p), one introduces an intertwiner to combine the eight 
spins associated with the four faces. Therefore, after integrating all 77's and £'s, the 
spins and the intertwiners combining the spins are left. The path integral is now 
written as a sum over spins and intertwiners representing surface-like excitations. 
The integrations of 77's, £'s and p's can be explicitly done by using the formulae 



and (flop. Integrate 77's and p's as follows. 



z As = fd( n n Evw(2v + vxuu^u), 



(31) 



xeA 3 0m/)=(i2), 1 

(20), (01) 



with 



A^(/3):= /'dp 
Jo 



(3-^1 dkkJ 2W+l (k)J 2lllu+1 (k) 



P~ 4 5(0) as /3->oo. 



(32) 



In the limit j3 — > 00, A Mi/ loses the dependence on and the path integral can 
be rewritten as follows. 



lim Z A3 = / d( J] II E(V + l)IXw(^) + Xw(-C^)], (33) 

(20), (01) 

up to an overall multiplicative constant. The first half is well known form for 3- 
dimensional lattice BF theory and is meant that spacetime is flat = 1. In ad- 
dition, we have another term also meaning that spacetime is flat but = — 1. 
The presence of the additional term is due to the definition of the lattice action. 
Because of the definition of the lattice action, the path integral contains a form anal- 
ogous to J dxe lxsind = 5(8) + 5(9 — n) instead of the well known form analogous to 
/ dxe %xd = 5(9). We have understood that the path integral of the model in 3 di- 
mensions can be reduced to (a generalization of) that of 3 dimensional lattice BF 
theory. 

4.2 Expectation values 

We compute the expectation values of four quantities, (r/ M ), (TrU^) and (Sa 3 ), 
and find all of them vanish. However, their meanings are different as discussed below. 
Remember that the path integral contains the integrations of 18 degrees of freedom 
per lattice point (9 for and 9 for r/ M and p M ). If one translates them into the 
canonical theory, 3 of and 3 of rj^ and p M are Lagrange multipliers imposing the 6 
constraints per lattice point. The other degrees of freedom are 6 degrees of freedom 
constrained and 6 on the constraint surface. The degrees of freedom constrained do 
not contribute to the path integral. When we compute the expectation values, the 
integrations over the degrees of freedom on the constraint surface make non trivial 
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contributions. All of the local degrees of freedom on the constraint surface are gauge 
degrees of freedom in the 3 dimensional case. 

(rjfj) and are basic variables of the theory. They are SU(2) gauge dependent 
quantities. It is known that the expectation value of basic variable of lattice QCD 
vanishes unless an appropriate gauge is fixed because there is a symmetric structure 
along the SU(2) gauge trajectories. The same Gauss gauge constraint structure is 
also present in general relativity in the present formulation. Therefore, (r/ M ) and 
should vanish if the model is consistent. 

Let us compute the expectation values of (77^) and Insert 770(2) into (|29f) and 
integrate T]q(x). This results in Z 2 i = I12 ±| associated with fn(x). On the other hand, 
the integrations of T]2(x), t]q(x — e 2 ) and 772 (a: — e°) result in l w = Z01 associated with 
foi{x), hi = I12 associated with fvz{x — e 2 ) and Z10 = Z01 associated with /oi(x — e°) 
respectively, fwix), foi(x), fn(x — £ 2 ) and foi(x — e°) share the edge where Ci(^) 
is defined. Hence, the integration of (i(x) introduces an intertwiner combining the 
spins associated with the four faces. However, no intertwiner consistent with these 
spins exists because the sum of the 8 spins is a half integer. Therefore, we conclude 
(Vo) — 0. In the same way, we find (771) = (772) = 0. 

In order to compute (Co), insert (o( x ) into ( |3T1) and integrate Co(^)- The edge 
where (o(x) is defined is shared by the four faces foi(x), fo2(x), foi(x — e 1 ) and 
fo2(% — £ 2 )- Since each of the four faces is associated with two equal spins after the 
integrations of 77's, the sum of the 8 spins associated with the four faces is an integer. 
There is no intertwiner combining an integer and | to produce the spin-0. Therefore, 
we conclude = 0. In the same way, we find (Ci) = (C2) = 0. 

(TrU^) is the expectation value of the Wilson loop and is an SU(2) gauge inde- 
pendent quantity. It is a good physical observable for Gauss gauge invariant theories 
such as lattice QCD. However, the Gauss gauge invariance is not enough to be a phys- 
ical observable for quantum gravity. Physical observables in quantum gravity must 
be invariant not only under SU(2) gauge transformations but also under spacetime 
diffeomorphisms. The Wilson loop is not invariant under diffeomorphisms. Therefore, 
it is possible that (TrU^) vanishes because of a symmetric structure along diffeomor- 
phism gauge trajectories unless an appropriate gauge is fixed. If it is the case, the 
vanishing result can be understood as a consequence of the fact that quantum gravity 
is a spacetime diffeomorphism invariant theory. 

The reasoning for (£ /t ) to vanish also holds for (TyU^) to vanish. Since every face 
is associated with two equal spins after the integrations of 77's, there is no way of 
combining the 8 spins associated with the four faces and spin | of TtU^ to produce 
the spin-0. Therefore, we find (TrU^) = 0. The fact that every face is associated 
with two spins is a notable difference of the model from lattice QCD. One of the 
reasons for this fact is the presence of the variable in addition to the connection 
variable A 1 ^ in the action. Note that (TrU^) = does not contradict to the fact that 
the same quantity evaluated with flat connection is not zero but one since this loop 
can be contracted to a point. This fact is restored by the computation of the quantity 
with a gauge fixing corresponding to either U^ v = 1 or — 1. The vanishing result is due 
to the symmetric structure consisting of U^ v = 1 and —1 on the constraint surface. 
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(Sa 3 ) is the expectation value of the action. It is invariant under SU(2) gauge 
transformations and spacetime diffeomorphisms. The value of the action on the con- 
straint surface is identically zero. Therefore, (Sa 3 ) must vanish if the path integral 
successfully eliminates the contributions from the degrees of freedom constrained. 
This must be the case if the path integral adopts the exponential oscillation form 
with f3 — > oo. On the other hand, if the path integral adopts the exponential decay 
form, then the dominant contribution comes from the value of the action far away 
from zero and hence (Sa 3 ) unlikely vanishes. This fact can be checked by replacing 
(3 in the path integral by i/3. 

In order to compute (Sa 3 ), take a derivative of ( |31| ) with respect to (3 (with i 



multiplied) instead of inserting S& 3 into (|29[) . This procedure results in the following 
quantity inserted in the path integral. 

E E ''^^ (34) 

(20), (01) 

Since there is no particular place or direction in purely empty spacetime, the every 
term must have equal contrubution and we compute one of the terms as follows. 

ij-^M _ At , 4*[J 2W+1 (2/3)] 2 _ U(, 1 _ 2roo tt, o; , ^ vr 



(35) 

as (3 goes to infinity. Here, we have used the asymptotic formula J m {x) — > cos(x — 
^ — |) for large x and the fact that A^ IU goes to /5 — 4 t5(0) or more precisely c/?~ 3 
as /3 — > 00, where c is a constant. We do not need the detailed value of c. This 
fact can be understood by rescaling (3 in (^2]). From (p5|) , we conclude (^3) — >■ as 
/5 — > 00. Note that if one replaces (3 by i(3 in the path integral, then one finds (j35D 
with the cosine replaced by the hyperbolic cosine. In this case, (Sa 3 ) diverges instead 
of converging to zero as (3 — >• 00. 



5 Conclusion 

In the present work, we constructed a surface-theoretic lattice quantum gravity model 
in 4 dimensional Riemannian spacetime based on the SU(2) Ashtekar formulation of 
general relativity. We introduced a version of the action of general relativity and 
defined its lattice version on a fixed hyper cubic lattice. We introduced a dimensionless 
"(inverse) coupling" constant so that the magnitude of the action is finite per lattice 
point. We defined a path integral whose integrand has the exponential oscillation 
form so that it eliminates the degrees of freedom constrained. The finiteness of the 
action with finite lattice allowed the path integral without gauge fixing finite. We 
expanded the path integral in the SU(2) characters and showed that the path integral 
can be written as a sum over surface-like excitations in spacetime. We showed that 
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the 3 dimensional version of the model exists and its path integral is reduced to 
that of 3-dimensional lattice BF theory. Therefore, we considered the model as a 4 
dimensional generalization of the Ponzano-Regge model with local degrees of freedom. 
We examined the expectation values of two basic variables of the theory, the Wilson 
loop and the action of the model in 3 dimensions and showed all of them vanish. We 
discussed the meaning of each of them and understood that the model has a chance 
of capturing the physical degrees of freedome of general relativity. 



A Useful formulae 

In the present work, we have used the following mathematical formulae. 

= £ 2 — ^■WaOttW, (36) 

3 X 

J m ( x ) = — T d6e ixsiae - ime . (37) 

The first formula can be easily proved and the second is a definition of the Bessel 
function of the first kind. Xj{U) is the character of SU(2) element U in the spin-j 
representation and j is a spin taking values 0, |, 1, | • • •■ 



(-l) n - m Diil(U)=D U l:_ m (U), (38) 



dUD^ n {U)D^* n ,{U) = (40) 

D^ n (U) is the spin-j representation matrix of SU(2) element U and m and n run from 

— j through ? with the increment 1. ( ^ 1 ^ 2 ^ ) is the so-called 3i-coefficient and 

V mi m 2 m J 

defined by 

( 31 32 3 '):= t^Li^0i^i;j2m 2 k-m), (41) 

with the Clebsch-Gordan coefficient (jifrii, j 2 m2\jm) . The asterisks mean the com- 
plex conjugate and the sum of j is taken over \j x — J2I through j\ + j 2 and the sums 
of n and m over —j through j. 
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